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Chapter 7 Cross product 


We are getting ready to study integration in several variables. Until now we have been 
doing only differential calculus. One outcome of this study will be our ability to compute 
volumes of interesting regions of R”. As preparation for this we shall learn in this chapter 
how to compute volumes of parallelepipeds in R?. In this material there is a close connection 
with the vector product, which we now discuss. 

















A. Definition of the cross product 











We begin with a simple but interesting problem. Let x, y be given vectors in R®: x = 
(11, 2,73) and y = (y1, y2, y3). Assume that x and y are linearly independent; in other words, 
0, x, y determine a unique plane. We then want to determine a nonzero vector z which is 
orthogonal to this plane. That is, we want to solve the equations xe z = 0 and yez =0. We 
are certain in advance that z will be uniquely determined up to a scalar factor. The equations 
in terms of the coordinates of z are 





L142 + £o2%q +2323 = O, 


Y121 + Yoze + 323 = 


It is no surprise that we have two equations but three “unknowns,” as we know z is not going 
to be unique. Since x and y are linearly independent, the matrix 


XZ, LQ IX 

Y1 Y2 YB 
has row rank equal to 2, and thus also has column rank 2. Thus it has two linearly independent 
columns. To be definite, suppose that the first two columns are independent; in other words, 


L1Y2 — L2Y1 F# 0. 


(This is all a special case of the general discussion in Section 6B.) Then we can solve the 
following two equations for the “unknowns” z, and 29: 


U1Z4+%9% = —1323, 


Y121 + Y2% —Y323- 


The result is 
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— 243 — ©3Y2 

Qi. tS Yeo se ea 
L1Y2 — L2Y1 
L3Y1 — L1Y3 
24.0 = 
L1Y2 — L2Y1 


Notice of course we have an undetermined scalar factor 23. 
Now we simply make the choice z3 = x1y2 — Xoy1(4 0). Then the vector z can be written 


4, = L2Y3 — @3Y2, 
22 = X3Y1 — L193, 
£30 = = X1Y2 — 21. 


This is precisely what we were trying to find, and we now simply make this a definition: 





DEFINITION. The cross product (or vector product) of two vectors x, y in R® is the vector 


LEXY = (Loy3 — L3Y2, T3Y1 — 13, Liy2 — Tet). 


DISCUSSION. 1. Our development was based on the assumption that x and y are linearly 
independent. But the definition still holds in the case of linear dependence, and produces 
x x y=0. Thus we can say immediately that 


x and y are linearly dependent <>2zxy = 0. 


2. We also made the working assumption that x,y2 — rey, # 0. Either of the other two 
choices of independent columns produces the same sort of result. This is clearly seen in the 
nice symmetry of the definition. 


3. The definition is actually quite easily memorized. Just realize that the first component of 
z=2X y is 

21 = ©2Y3 — ©3Ya2 
and then a cyclic permutation 1 — 2 — 3 — 1 of the indices automatically produces the other 
two components. 


4. A convenient mnemonic for the definition is given by the formal determinant expression, 


1 fj &k 
xxy = det| x7, 2 23 
4. Y2 YB 
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In this expression the entries in the first row are the standard unit coordinate vectors, and 
the “determinant” is to be calculated by expansion by the minors along the first row. 


5. We are left with one true ambiguity in the definition, and that is which sign to take. In our 
development we chose 23 = ©1Y2 — ®2y1, but we could have of course chosen z3 = Loy — 11Ye. 


In this case, the entire mathematical community agrees with the choice we have made. 


6. Nice special cases: 


ixj = k, 
jxk = 
kxt = j. 





7. By the way, two vectors in R® have a dot product (a scalar) and a cross product (a 
vector). The words “dot” and “cross” are somehow weaker than “scalar” and “vector,” but 


they have stuck. 
ALGEBRAIC PROPERTIES. The cross product is linear in each factor, so we have for 


example for vectors x, y, U, UV, 
(az + by) x (cut+ dv) = acz x u+adz x v+ bey x ut bdy x v. 


It is anticommutative: 


yxXu = -£X y. 
It is not associative: for instance, 
ix(ixj) = txk = -j, 
(ixi)xj = 0x7 = 0. 











PROBLEM 7-1. Let x € R® be thought of as fixed. Then x x y is a linear function 
from R* to R® and thus can be represented in a unique way as a matrix times the column 
vector y. Show that in fact 























ex y= | 2% O. =f | ws 
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PROBLEM 7-2. = Assuming x ¥ 0 in the preceding problem, find the characteristic 
polynomial of the 3 x 3 matrix given there. What are its eigenvalues? 


B. The norm of the cross product 


The approach I want to take here goes back to the Schwarz inequality on p. 1-15, for 
which we are now going to give an entirely different proof. Suppose then that x, y € IR". We 
are going to prove |xe y| < ||z|| ||y|| by calculating the difference of the squares of the two 
sides, as follows: 














nm n n 2 
lll? vi? -(wey)?? = Soa? Soy? - (Soom) 
i=l i=l i=1 
= oy = Saini >) ayy; 
i=l j=l i=l = 
= Ds, iy; _ » LEYiX jYj 


ij=l j=l 
= Dosing + Do aiut + 0 aia; 
i<j i=1 i>j 
_ S° LzYiL 5 Yj — oy = y LiYil jf Yj 
i<j i=] i>j 
= Sosiyg+ So apy} -25— siyayy; 
i<j j>i t<j 
= So (afyf — 2riysnyyj + 25y?) 
i<j 
= > (xiyy — iii) 
i<j 


This then proves that ||zx||? ||y||? — (ey)? > 0. 
If we specialize to R®, we have 





lz II? lly Il? —(xe y) = (1y2 — reyi)” + (x2y3 — r3Yy2)” + (z1y3 — T3y1)°. 


But the right side is precisely 23 + 27 + 23 in the above notation. Thus we have proved the 
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wonderful Lagrange’s identity 





lil Ill? = @ey)? + lle x yl. 











This identity relates norms, dot products, and cross products. 


In terms of the angle 6 between x and y, we have from p. 1-17 the formula 7 ey = 
||| ||y|| cos 6. Thus, 








lz x yll = [ell llyllsin 8. 








This result completes the geometric description of the cross product, up to sign. The 
vector x x y is orthogonal to the plane determined by 0, x and y, and its norm is given by the 
formula we have just derived. An even more geometric way of saying this is that the norm of 
x x y is the area of the parallelogram with vertices 0, x, y, x + y: 


@Xxy 


xXt+y 
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PROBLEM 7-3. The Lagrange identity can easily be generalized. Using the following 
outline, prove that for x, y, u, v € R®, 














(xx y)e(uxv) = (reu)(yeu) —(rev)(yeu). 


Outline: the right side equals 


(Eee) (Es) - (Eas) (Toe 


= > Ly UU; = 55) (why?) 
1) 

= 2 Li Yi(UzVi = UjV;) (why?) 
a) 
| 

= 5 Yo (inj — 2jy:)(wivj — ujvi) (why?) 

ij 

= So (xiyj — 25y:)(uivy — ujvs) (why?) 

<j 


= (¢xy) elu xv). 


C. Triple products 


The formal representation of x x y as a determinant enables us to calculate inner products 
readily: for any u € R? 





+ 9 & 
(cxy)eu = det {x 13) ° (wi + Ugj + ugh) 
Y. Y2 YB 
Uy U2 UZ 
= det|{z, wr 23 
M1 Y2 Ys 


(Notice how clearly this formula displays the fact that (a x y)ex = (a x y)ey =0.) 
This number (x x y) eu is called the scalar triple product of the vectors x, y, and u. The 
formula we have just given makes it clear that 
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(xcxy)eu = (uxaujey = (yxujer 
= -(yxz)eu = -(xxujey = —(uxy)eg. 
Thus we have the interesting phenomenon that writing x, y, u in order gives 
(cx y)eu = re(yxu). 
So the triple product doesn’t care which you call x and which you call e. For this reason, it’s 
frequently written 
[x,y,u] = (4x y)eu. 


PROBLEM 7-4. The vector triple product is (x x y) x u. It can be related to dot 
products by the identity 


(x xy) X u=(reu)y— (yeu)z. 


Prove this by using Problem 7-3 to calculate the dot product of each side of the proposed 
formula with an arbitrary v € R°. 














PROBLEM 7-5. Prove quickly that the other vector triple product satisfies 


ax(yxu) = (weu)y—(vey)u. 


The identities in Problems 7—4 and 7-5 can be remembered by expressing the right side of 


each as 
vector triple product = (outerefar) near — (outerenear) far. 


The scalar triple product is exactly the device that enables us to compute volumes of 
R°. Consider three vectors x, y, u in R®. They generate a region D called 








parallelepipeds in | 
a parallelepiped as follows: 


D = {oa by on |< a <1, 05541, 0 <e6= 1}, 





PROBLEM 7-6. Prove that 


[cx y,y x u,u X 2] = [2,y, ul’. 
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We use the usual intuitive definition of volume as the area of the base times the altitude. If we 
regard x and y as spanning the parallelogram which is the base, its area is of course ||x x y||. 
The altitude is then calculated as the norm of the projection of the side u onto the direction 
orthogonal to the base; of course, x x y has the required direction: 


XXy 





plane of 0, x, y 





In terms of the angle 0 between x x y and u, the altitude is 
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x 
ul] [cos8] = full HEX woul 
le x all Tel 
. Nesaien 
le x4l 


Thus the volume of D is 
vol(D) = |(a xy) eul. 
It is therefore computed as the absolute value of the scalar triple product of the three spanning 


edges of the parallelepiped. 
In terms of the determinant, 


vol(D) = | det XZ LQ XM 


Yo Yo YB 





U; Ug U3 


Notice in particular that vol(D) = 0 <=> the vectors z, y, u are linearly dependent <=> the 
parallelepiped D lies in a two dimensional subspace of R”. 














D. Orientation 


We now have the cross product characterized geometrically, up to a sign. Given linearly 
independent x, y € R®, the cross product has direction orthogonal to both x and y, and norm 
equal to ||z|| ||y||sin@, where @ is the angle between the vectors x and y. There are exactly 
two vectors with these properties, and we now want to discuss a geometric way to single out 
the cross product. 

The concept we need is important in general dimensions, so we shall deal with R” for any 








n. 


DEFINITION. A frame for R” is an ordered basis. That is, we have n linearly independent 
vectors (1,---,; Qn, and they are written in a particular order. For example, 


1059; ky, {@ a 2), 4; —k}, {2, kj} 


are three distinct frames for R®. 











DEFINITION. Let {y1,...,y,} and {v,...,%,} be frames for R”. Then there exists a 
unique n x n matrix A such that 
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This matrix is necessarily invertible. We say the two frames have the same orientation if 
det A > 0; and opposite orientation if det A < 0. 


PROBLEM 7-7. In the above definition we have used the column vector representa- 
tion of vectors. We may then write the corresponding n x n matrices 


@ = ($1---Yn), 
a (Bis. Py). 


Prove that the frames have the same orientation — > det ® and det WV have the same 
sign. 


PROBLEM 7-8. Prove that “having the same orientation” is an equivalence relation 
for frames. (Thus, as far as orientation is concerned, there are just two types of frames.) 


PROBLEM 7-9. Prove that if {y1,..., gn} and {v1,...,Wn} happen to be orthonor- 
mal frames, the matrix A in the definition is in O(n); and the frames have the same 
orientation = > A € SO(n). 


Let us say that the coordinate frame {é,,...,é,} (also written {7,j,k} in the case of 
R*) has the standard orientation. The corresponding matrix is the identity matrix J, and 
det J = 1. Thus we say that an arbitrary frame {y1,...,¢,} has the standard orientation if 
the corresponding matrix ® = (Y~1... Yn) has det ® > 0. 





PROBLEM 7-10. For which dimensions do both {¥1, ¢2,..-, Yn} 
and {Yn, ¥1, Y2,---;Yn-1} have the same orientation? 


PROBLEM 7-11. For which dimensions do {y, Yo,.-., Qn} and 
its “reversal” {yp,..., 2,1} have the same orientation? 














Now we consider R? and two vectors x, y whose cross product «xy # 0. Writing everything 
as column vectors, we then have from Section C that 


det(x ya xy) = |lxxyll? > 0. 


Thus the frame {z,y, x x y} has the standard orientation. This fact is what settles the choice 
of direction of the cross product. 
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SUMMARY. The cross product x x y is uniquely determined by the geometric description: 
(1) x x y is orthogonal to x and y, 
(2) ||a x y|| = area of the parallelogram with edges x and y, 
(3) in case x x y £0, the frame {z, y, x x y} has the standard orientation. 


Since the cross product is characterized geometrically, we would expect the action of ele- 
ments of O(3) to have a predictable outcome. Indeed this is so: 


THEOREM. Let A € 0(3). Then 
A(z x y) = (det A)Ax x Ay. 


Thus in case A € SO(3), then A(x x y) = Ax x Ay. 


PROOF. We use the characterization in terms of column vectors 





(cx y)eu = det(zx y u), all u € R®. 


Then for any 3 x 3 real matrix A, 


(Ar x Ay)e Au = det(Ax Ay Au) 
= det A(z y u) 
= det Adet(z y u) 
= det A(x x y) eu. 


But 
(Ax x Ay)e Au = A‘(Ax x Ay) eu. 





Since we have two vectors with the same inner product with all u € R®, we thus conclude 
A‘(Ax x Ay) = det A(x x y). 
If A is also invertible, then 
Az x Ay = det A(A*)"'(z x y). 


Finally, if A € O(3), then (A*)~' = A and det A = +1. 
QED 
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The outcome of our discussion is that the cross product enjoys the double-edged properties 
of being characterized completely geometrically and of being easily computed in coordinates. 
In Chapter 1 we experienced the same happy situation with dot products on R". The dot 
product x e y can be completely characterized as |||! ||y||cos@, where 6 is the angle between 
x and y, but also it is easily computed in coordinates. 














3 


xX 


xy 


E. Right hand rule 


This section just gives another way of presenting orientation. Given a Euclidean three- 
dimensional space, we may think of first choosing an origin and then placing an orthogonal 
coordinate system 271, %2, 73. We say this is a right-handed coordinate system if when you 
place your right hand so that your fingers curl from the positive x-axis through 90 degrees 
to the positive x2-axis, your thumb points in the direction of the positive x3-axis. 


Another way of saying this is to place your right hand so that your fingers curl from 7 to 
j. Then your thumb points in the direction k if the coordinate system is right-handed. 


Notice that 7 x j = k. This then gives the right-hand rule for cross products, supposing 
we have a right-handed coordinate system. Supposing that x x y 4 0, place your right hand 
so that your fingers curl from x toward y through an acute angle. Then your thumb points in 
the direction of x x y. 


This is easily proved by the device of using SO(3) to move your right hand from the {?, j, i} 
frame to an orthonormal frame associated with the frame {z,y,2 x y}. Here’s the picture, 
looking “down” at the plane determined by 0, 2, y: 
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y’ 
@ @ 
oO x’ x 
Let 2’ = x/|\x||, and let y’ be the unit vector shown: 
! y—aeya' 
i= / i 
lly—aleya'| 
Thus y = c2’ + coy’, where cy > 0. Then 
txy =] Gexy = elle a xy. 


Thus x x y is a positive scalar multiple of x’ x y’. Now define the matrix A in terms of its 
columns, 

A= (ey a x4): 
Then Ai = 2’, Aj = y', Ak = x’ x y’. And we know A € SO(3). Thus A takes the right- 
handed frame {2, j, k} to the right-handed frame {2’, y’, 2’ x y’}. This concludes the proof of 
the right-hand rule. 


F. What is SO(3)? 


This seems to be a good place to discuss the fine details of SO(3), the special orthogonal 
group on R°. The members of this group are of course 3 x 3 orthogonal matrices with determi- 
nants equal to 1. Let A € SO(3). This matrix of course “acts on” R? by matrix multiplication, 
sending column vectors x to the column vectors Ax. We are going to prove that this “action” 
can be described as an ordinary rotation of R? about an axis. 

We first pause to describe rotations. Suppose that w is a fixed unit vector in R?, and that 
6 is a fixed angle. We use the standard right-handed orientation for R?. Now we take a view 
of R? looking from w towards 0, and rotate the entire space counter-clockwise about the axis 
w through the angle 6. Let us call the resulting matrix R(w,@). A couple of things should 
be noticed. One is that a counter-clockwise rotation with a negative angle @ is actually a 
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clockwise rotation with the positive angle —@. The other is that adding any integer multiple 
of 27 to @ results in the same rotation, according to our convention here. 

Now we show how to calculate R(w,6@). To do this use a right-handed orthonormal frame 
{t, 6,0}. This can be computed by choosing any unit vector t which is orthogonal to w, and 
then 0 is uniquely determined by the required formula + = w x u. Now think of projecting 
R® orthogonally onto the @, 6 plane, and then rotating through the angle 6. Here’s the view 
looking down from w: 











p ° v 
a 
Q e 
R(W, 6) p 
Thus we have 
R(w,0)u = ticosé+idsin8, 
R(w,0)o = —tsind+tcosé. 
Of course, since wW is the axis of rotation, 
R(w,O)w = w. 


The summary of this information is expressed in terms of matrices as 
R(w,O)(@bw) = (tcosO+tsind —tsind+tcosd w) 
cos? —sin@ 0 
(ad w) | sind cosd 0 
0 0 1 


Therefore 
cos@ —sind 0 


R(w,0) = (@0w){sin@ cosd 0] (@0 Ww). (x) 
0 O 4 
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Remember: the computation of the inverse of the orthogonal matrix (t 6 w) is just a 
matter of writing down its transpose 


ity Ss 


& 


Based on our understanding of what R(w,6@) is, we are certain that this formula does not 
really depend on the choice of &% and t. See Problem 7-16 for the explicit verification of this 
fact. 


You should notice that the inner matrix in the above formula is precisely R(k, 0). So the 
formula displays how we are able to move from a rotation about the x3-axis to a rotation 
about the axis determined by w. 


PROBLEM 7-12. Prove that R(w,6) € SO(8). 


PROBLEM 7-13. Using the above formula, show directly that 
Riw,0) =. 1: 


Then give a heuristic explanation of this equation. 


PROBLEM 7-14. Using the above formula, show directly that 


Then give a heuristic explanation of this equation. Also explain why 
we may thus assume without loss of generality that 0 <0 <7. 


PROBLEM 7-15. Using the above formula, calculate directly that 
R(w,0,)R(w, 02) = R(w, 6, + 2). 


Then give a heuristic explanation of this equation. 
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PROBLEM 7-16. Calculate the product and thus obtain from (x) 
0 —W3 W92 
R(w,0) = (cosé)I + (1—cos0)(ww;)+sind| ws OO uw 


—W2 Wy 0 


(This displays R(w,0) in a form that shows its independence of the 
choice of t and @.) 


PROBLEM 7-17. Show that the action of R(w,0) can also be expressed in the form 


R(w,0)x = (cosé)x + (1 — cos6)(we x)w +t (sin d)w x a, 


for all x € R®. Do this two ways: 
a. by simply using Problem 7-16; 


b. by describing the action of R(w, 0) geometrically. 


PROBLEM 7-18. Prove that for all xe w = 0, 


Give a geometric interpretation of this result. Conversely, suppose 
that for a given rotation there exists a nonzero x € R® such that 


Prove that 9 = 7 and xew = 0. (Of course, 6 could also be 7+ 2ma 
for any integer m.) 





PROBLEM 7-19. Prove that R(w,@) is a symmetric matrix <> 6 = 0 or 7. 
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PROBLEM 7-20. Backward Andy all his life has insisted upon using a left-handed 
coordinate system for R*?. What formula must he use for R(w,@) that is the analog of 
ours found in Problem 7-16? 














We now turn to the perhaps surprising converse of the above material, which states that 
the only matrices in SO(3) are rotations! The crucial background for this result is contained 
in the following 


THEOREM. Let A € SO(n), where n is odd. Then 1 is an eigenvalue of A. 
PROOF. Since A7! = A‘, we have 


det(A—I) = det(A— AA‘) 
= det(A(I — A’)) 
= det Adet(I — A‘) 
= det(I — A’) 
= det(I — A) 
= (—-1)"det(A-— J). 


Since n is odd, det(A — J) =0. 
QED 


PROBLEM 7-21. — Give a different proof of this result for the case n = 3, based on 
the characteristic polynomial 


det(A — AL) = —A2+ ad + BA+1. 





PROBLEM 7-22. Let A € O(n), but A ¢ SO(n). Prove that —1 is an eigenvalue of 
A. 


THEOREM. Let A € SO(3). Then there exists a unit vector w € R® and 0 < 6 < 7 such 
that A = R(w, 6). 





PROOF. Since A € SO(3), we know from the preceding theorem that 1 is an eigenvalue of 
A. Therefore, there exists a unit vector w such that Aw = w. 
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Choose an orthonormal basis &, 6 of the subspace of R® which is orthogonal to w; arrange 
them so that the frame {t,0,w} has the standard orientation: 


yw = UX v, 
oe = wx 4, 
U = UXW. 
Then 
Atiew = Atte Aw 
= uew (A is in O(3)) 
= 0 


? 


so At is also orthogonal to w. Thus Aw is a linear combination of t and #6: 


Ati = cy + C20. 


Since ||Aé|| = ||a|| = 1, c? +c$ = 1. Thus there exists an angle @ such that c, = cos0, 
co = sind. That is, 


At =tcosé+ sind. 


We conclude that 


At = A(wx&) 
= Aw x At (Section D) 
= wx (icosé+ dsin 8) 


= bcosd—itsing. 


We now see that A and the rotation R(w, 6) produce the same result on the frame {t, 6, w}! 
Since every x € R? is a linear combination of &, 6, w, we conclude that 








Ar = R(w,0)x for all x € R®. 


Thus A = R(w, 8). 
The final item, that we may require 0 < 0 < 7, follows immediately from Problem 7-14. 
QED 
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DISCUSSION. For a given A € SO(3), the representation A = R(w,6), where wW is a unit 
vector and 0 < @ < 7, is virtually unique. To see this we start with the result of Problem 7-16, 
and compute the trace of each matrix: 





traceA = 3cosé+(1—cos6)(w7 + ws + w3) 
= 3cos@+ (1—cos6@) 
= 2cos?+1. 
Thus 
traceA — 1 
cos? = —— 7 = 


Therefore, as 0 < 6 < 7, the angle is completely determined by the knowledge of traceA. 
Next notice that Problem 7—16 implies 


0 —W3 Wa 
Am A = 2sin0 W3 0 UW 
—W2 W1 0 


We conclude immediately that if 0 < 6 < a, the vector w is also uniquely determined. 
The only cases of nonuniqueness are thus when sin@ = 0, or, equivalently, A = A‘. If 
6 =0, then A = R(w,0) = J and the “axis” w is irrelevant. If 9 = 7, then we have 


A => =) 2). 


It is easy to check in this case that w is uniquely determined up to sign. This is because w 4 0 
and so some coordinate w; 4 0 and we obtain from the above equation 


2w; = j column of (A+ J). 


This determines w up to a scalar factor, and the fact that ||2|| = 1 shows that the only 
possibilities are w and —w. 

We now have a truly intuitive understanding of SO(3). In particular, the right-hand rule 
makes wonderful sense in this context, as we now explain. Suppose {%1, G2, 63} is an arbitrary 
orthonormal frame with standard orientation, and form the usual matrix ® with those columns, 


b = (1 Ge $3). 


Then © € SO(3) and ®' maps the frame {$1, G2, G3} to {i, j, 
R(w,0). Thus if you place your right hand to lie in the correct k 
your hand with R(w,6@) produces the ~3 = 1 X 2 position. 


k}. Moreover, ® is a rotation 
= 7X 7 position, then rotating 
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I would say it is almost amazing that a single rotation can move {(1, G2, 63} to {2, J, k}! 
Here is a numerical example. Let 


=] g 
A= ;]0 0 -l 
1 0 0 


Clearly, A € SO(3). Then we compute traceA = 0, so that cos@ = —} and we have @ = 27/3. 
Then 


0 —W3 W2 
A= = 2sin60 W3 0 Wi 
—W2 Wi 0 
becomes the equation 
0 -1 -l 0 —W3 Wa 
1 0 -1] = vV3{[ vw O —w], 
14a 0 —w, wi 0 
so we read off immediately that 
1 


Thus 


PROBLEM 7-23. Let 











1 1 1 
y Y Vie 
v3 v2 v6 
vi (8 


Find w and @ such that A = R(w, 0). 


PROBLEM 7-24. Compute 


R(d, 0)" = R(?,?). 
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PROBLEM 7-25. — Show that the given matrix is in SO(3) and calculate the unique 
w € R? and 0 < 0 < o for it: 














1 as =20 28 
mT 28 30 —-4 
—20 20 35 


G. Orientation reversal 

It is a relatively easy task to understand those matrices in O(n) which are not in SO(n). 
Suppose A is such a matrix. Of course, det A = —1. We conclude from Problem 7-22 that 
—1 is an eigenvalue of A. Therefore, there exists a unit vector w: 


Aw = —w. 


Any unit vector w can be used to define a “reflection” matrix S' in that direction. Here’s 
a schematic description: 








ex 
Ww ! 
0 : 
@ ox 
The algebra is governed by 
(Str+arjew = 0, 
Sx = x+tw forsome tER. 
Solving these equations gives t = —2x% e w so that 


St=x2—2A(rcow)w. 


22 Chapter 7 


This matrix S is called “reflection in the direction w,” and we denote it S(w). 


PROBLEM 7-26. Prove that the matrix S = S(w) satisfies 


a. SW = —-w, 
by 8? =, 

ec. S € O(n), 
d. det S = —1. 


PROBLEM 7-27. Prove that 





Now we return to the situation of a matrix A € O(n) with det A = —1. We know from 
Problem 7-22 that there is a unit vector w satisfying Aw = —w. Now just define B = S(w)A. 
Clearly, B € SO(n). Thus, 

A = S(W)B. 
Thus A is the product of a reflection S(#) in a direction w and a matrix B € SO(n) with 
Bw = w. 


PROBLEM 7-28. It made no difference whether we chose to write B = SA or 
B = AS, for SA = AS. Prove this by noting that SAw = ASw and then considering 
those vectors x which are orthogonal to w. 





In particular, for n = 3 the matrix A has the representation 


A = S(w)R(wW,0) = R(w,0)S(w). 


PROBLEM 7-29. Show that in this n = 3 situation 


0 —W3 


A = (cos@)I — (1+ cos@)(wjw;)+sin6| wy 0 —wi 
—W2 Wy, 0 
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PROBLEM 7-31. Refer back to inversion in the unit sphere in R”, p. 6-34: 


x 


f(x) = Wap: 
|x|? 
Prove that 


bf(e) =[lelr?s (727). 





Thus we say that inversion in the unit sphere is an orientation—reversing conformal mapping. 


PROBLEM 7-32. Let wu, v be linearly independent vectors in R? and let u’, v’ also be 
vectors in R?. 


a. Suppose that A € SO(3) satisfies 


Prove that u’ and v’ are linearly independent and 


al ale 


al llo"ll, (*) 


uev = uev’. 


b. Conversely, suppose that (*) is satisfied. Prove that u’ and v’ are linearly indepen- 
dent, and that there exists a unique A € SO(3) such that 


(HINT: Ala ou xe) = vo wx v),) 
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PROBLEM 7-33. Under the assumptions of the preceding problem prove that 


(uvuxv)t=luxv||A(v x (uxv) -—ux(uxv) ux)’. 





